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1.0 Introduction 

The ARC (Altimeter Residual Computation) computer 
program was developed at the Wallops Flight Center, Wallops 
Island, Virignia,by Jim McMillan for Wolf Research and 
Development Group. The program acts as the final altimeter 
preprocessor before the data is reformatted for external 
users. It calculates all parameters necessary for the 
computation of the altimeter observation residuals and the 
sea surface height. 
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2.0 ARC Input Files 

The following pages contain a description of all of the 
ARC input files. The user supplied data control cards will 
be described in Section 4. 

2.1 GEODYN RV Tape 

The GEODYN RV tape is a binary, unformatted, 800 BPI, 
standard label tape which supplies ARC with the satellite 
latitude, longitude, and height above the reference ellipsoid 
at fixed tine intervals for interpolation. The tape format is 
given in Table 2-1. It should be noted that the only 
quantities on the tape which are used by ARC are DAY, SATLAT, 
SATLON, and SATHGT. ARC reads the GEODYN RV tape from file 
code 1. 

IMPORTANT: The ARC program is designed to interpolate 

the RV tape and cannot perform extrapolation. 
Therefore, the RV tape must contain at 
least NORDER/2 records of data prior to the 
first observation and after the last 
observation to be processed, where NORDER 
is defined as the interpolation order (see 
Section 4.2.2). 

2.2 Geoid Permanent File 

The Marsh-Vincent geoid is supplied to the ARC computer 
program by means of a random disk permanent file with a 


2.1 


constant record size of 133 words. The data is valid only 
for latitudes greater than -66° and less than +66°. If an 
attempt is made to evaluate the geoid outside of these limits, 
the message "SATELLITE LATITUDE - __IS OUTSIDE RANGE OF GEOID - 
STOP," and the program will terminate. The data has been 
written on the disk in a one-degree latitude by one-degree 
longitude grid. Each disk record consist' of a vector of 
all of the geoid values for a given longitude. ARC evaluates 
the geoid by linear interpolation at one-second intervals 
when computing observation residuals. The value of the geoid 
written on the output tape is evaluated at the time of the 
smoothed observation. ARC reads the geoid permanent file from 
file code 2. 

2.3 Data Tape 

The ARC input data tape, which is normally supplied as 
an output of the GAP computer program, is a binary, unformatted, 
800 BPI, standard label tape which supplies ARC with the altimeter 
data to be processed. The formats for the header record and the 
low and high data rate data records are given in Tables 2-2 to 
2-4. It should be noted that the only quantities on the tape 
which are used by ARC are MJDATE, FODAY, SMOALT, NWORDS, 

FRAMTI , DOY, YEAR, and CALT (or ALT). Arc reads the data tape 
from file code 29. 


2.2 


Word 

Variable Name 

Type 

Description 

1-2 

DAY 

D 

Days from Jan 0.0 of 
the reference year. 

3 

IYMD 

I 

Date in YYMMDD. 

4 

IHM 

I 

Hours and minutes in 
HHMM. 

5 

SEC 

R 

Seconds . 

6-7 

POSVEL 

D 

Inertial position and 
velocity. 

18-19 

SATLAT 

D 

Satellite latitude in 
degrees . 

20-21 

SATLON 

D 

Satellite longitude in 
degrees . 

22-23 

SATHGT 

D 

Satellite altitude in 
meters. 


TABLE 2.1. The 

GEODYN-RV 

Format 
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3.0 ARC Output Files 

The following pages contain a description of all of the 
ARC output files except scratch disk files and SYSOUT files. 

3.1 Plot Tape 

The plot tape generated by ARC is a non-standard label, 
556 BPI, "pen instruction" CalComp plot tape with a buffer 
size of 1024 words, suitable for mounting on the WFC CalComp 
763 plotter. The tape is not normally generated but can be 
created by the appropriate data control card ( see Section 
4.2.4). ARC writes the plot tape file code 3. 

IMPORTANT: Since ARC generates one frame of plot 

(16 inches) for every minute of data, 
it is usually not desirable to request 
that the tape be plotted until it is 
known which minutes are of interest. 
Therefore, the program outputs a minute 
by minute summary of the block numbers 
of the plot tape. The user can then 
request any segment of the tape which 
he needs to have plotted by merely 
noting the block numbers on the computer 
request form and checking PLOT. 

3.2 Output Tape 

The ARC output data tape, which is normally used as 
the input to the CALTOR computer program, is a binary. 


unformatted, 800 BPI, standard label tape. The formats for 
the header record and the low and high data rate data record 
are given in Tables 3-1 to 3-3. The parameters which ARC 
changes when writing the output tape are PREPIN, LAT, LON, 
ABIAS , TREF, GKITE, THITE, SATHT , CSIGH, SSHITE, and SMOSSH. 
ARC writes the output data tape on file code 4. 
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4.0 Data Control Cards 

The A?„C computer program has been designed to operate 
without user supplied data cards whenever possible. However 
there may be some applications which require special execution 
directives. These directives are supplied to the program by 
means of the optional data control cards. 



4.1 Mandatory Data Control Cards 

4.1.1 Title Card 

The first data control card encountered by the ARC 
program is considered as the run title and is printed on every 
page of output, on the run summary, on the printer plots, and 
on the CalComp plots. The title is alphanumeric, 80 characters 
long, and should be centered on the data card. 

4.2 Optional Data Control Cards 

All optional data control cards are read according 
to the following format: 


FORMAT (AS,I5,3F10.5,A40) 



where the first field is the card identifier, the next four 
fields supply executive information, and the last field is an 
unused field which may contain comments, numbers, etc. If an 
optional data control card with an unrecognizable identifier 
is encountered, the program will terminate. 

4.1 


4 . 2.1 *BIAS Control Card 

The *BIAS control card is used to supply the program 
with the value of the altimeter bias. If no card is present, 
the bias will re set to zero. 


•♦•BIAS 

III 


OODDi OBODOOOOOOCDOCeOJSOCCOOOGnOODCi'OOGOOOOCOOOGOOOOOBOOOOODOEDaOOGOOCOGOODGCOOO 

I S 3 4 5 9 > ft f 13 If !2 13 K ISM IMS li 25 3 2?»34 25 25 27 112133 21 32 3135 35 33 3/ 3J33 <3 41 <2 43 44 45 4*43 <3 4*35 5? 57 !3 si Si Si JJ *3 53 53 S! £! S3 £t 1-5-5 5? 63 S3 70 71 72 7* 25 7* 71 71 13 84 

ml i n i m 1 1 1 m i ; i n i 1 1 1 1| n 1 1 n 1 1 1 m i n 1 1 n n i m n 1 1 1 1 ii 1 1 1 1 n n 1 1 1 1 m n 1 1 
2 | 2 21 22 2 2 2 2 222 2 2 2 2 2 | 2 2 2 2 2 2 2 2 21 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 
3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 31 3 3 3 3 3 3 3 3 3| 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
1-444444444444444 444444<444<44<4<4 4 <<4444<4<44444444444<44C4444444<444A4444444444 
5 5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5| 5555555555555555555555555555555 5 5555555555 5 5555555355555555555 
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I 2 5 4 S C 7 3 • i9 t: 1! ft U ’.5 1> W •: -- l 22 2; :: li 2 1 It 2i 20 »i 5? 27 24 Si5S 27 :i 3) <3 x < <2 43 44 « V 13 5* 52 52 :* 53 *i 51 aSCi Si *2 EJ i< C. e$ 57 f 2 iZ 73 7i Ji jj ?1 Ji 7i ; ; li #9 79 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*BIAS 

1-5 

A5 


BIAS1 

11-20 

F10 . 5 

global mode bias 

BIAS2 

21-30 

F10 . 5 

intensive mode bias 



4.2 


4.2.2 “EDIT Control Card 

The *EDIT control card is used to supply the program 
with the -editting criterion for the plotting of the altimeter 
residuals. It is necessary to edit out bad residuals since 
ARC scales its plots based upon the maximum and minimum values 
to be plotted. If no card is present, observation residuals 
greater than 50 meters in absolute value will be editted. 
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1 000 . 0 
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NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*EDIT 

1-5 

A5 


RED IT 

11-20 

F10.5 

observation residual 


editting criterion 


4.3 


4.2.3 *ORDR Control Card 

The *0RDR control card supplies the program with the RV tape inter- 
polation order. The maximum acceptable value is 7 and the minimum accept- 
able value is 3. It should be noted that the interpolation order effects 
the number of records needed on the RV tape before the first observation 
and after the last observation (see Section 2,1), The order of the 
interpolation chosen should reflect the mesh size of the RV tape. If no 
card is present, the interpolation order will be set to 7. As of 2 Sept- 
ember 1976, the spacing on RV tapes was 30 seconds for tapes generated at 
WFC and 60 seconds for tapes generated at NWL. 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*0RDR 

1-5 

A5 


NORDER 

6-10 

15 

RV interpolation order 


4.4 


2 September 1976 


4.2.4 *PLOT Control Card 


The *PL07 control card supplies the program with the 


plot option chosen by the user. If no card is present, the 
program will produce only printer plots for every minute of 
data. It should be noted that an output tape can be generated 

by the program when the user supplies a *PLOT control card. ‘ 

-PLOT 3 
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NAME 

COLUMNS 

FORMAT 

*PLOT 

1-5 

A5 

IPLT 

6-10 

15 
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DESCRIPTION 


$&AI Pa 


plot option 

0 = no plots 

1 = printer plots 

2 = printer and 

CalComp plots 

3 = printer and 

CalComp summary 
plots only 


4.5 


4.2.5 -Si'C? Control Card 


The ft ST0P control card is used to supply the program with 
the desired stop time. If no card is present, the program will 
terminate upon reaching an end of file on the observation tape. 

*STOP ' 740531 , C 115534.0 \ 
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NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*STOP 

1-5 

AS 


YMI) 

11-20 

FI 0.5 

YYMMDD. of stop time 

HMS 

21-30 

F10.5 

HHMMSS.SSS of stop time 
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4.2.6 *STR: Control Card 


The *STRT control card is used to supply the program with 
the desired start time. If no card is present, the. program will 
start with the first observation on the data tape. 


•*'STRT 74053 1.0 11 5534 . 0 

I I 


D|| 0| OCOPOOOOOl ccoo; OCSilDODCOi 0 0 C 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 41 0 0 0 D 3 0 0 0 0 0 0 3 0 0 3 0 
i i ] < l • ; i i io ii i: i) v. :i r n » up :: a mim 2 i»: 9 11 aim miKinuju 53 it sisistsisisr s» n*a » 11 a n n nr, :i nnu 

11 n 1 1 11 m n 1 ml 1 1 i :M n 1 m 1 1 m 1 1 it 1 1 1 m n 11 1 tt 1 1 1 1 1 m 1 1 1 1 1 1 1 1 11 1 11 m 1 1 1 1 

2 ! 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 'I 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

3 3j 3| 3 3 3 3 3 3 3 3 3 3 31 3| 3 3 3 3 3 3 31 3| 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

| 4 <44<44444< 4| < 4 4 4 4 4 4 4 4 4 4 4 4| 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 < 4 4 4 4 < 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 

5 5 5 5 5 5 5 5 5 5 5 5 5 5 5| 5 5 5 5 5 5 5 .i| | 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 a 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 

. e 5 C 6 6 C P E 6 5 6 5 £ 6 5 6 £ S 8 6 C 5 S £ S 5 S S c 6 C 6 E 6 6 £ C 5 6 6 6 6 C t 6 6 C- 5 E 6 C 6 C S E 0 S 6 6 5 t R 6 £ S E 6 £ 3 8 S 6 E 6 6 6 E £ S 6 

' 7 7 7 7 7 7 ) 7 7 7 7 ?| 7 7 7 3 7 7 7 7 7 7 7 ? 7 1 7 7 7 7 7 7 7 7 7 J 7 / 7 7 7 7 7 7 / 7 7 7 7 7 7 7 7 7 J 7 7 7 7 7 7 7 ? 3 7 7 7 7 1 1 7 7 7 1 7 7 1 7 7 

| £ t £ 8 B 5 8 g 8 8 8 E B S 8 E S| 0 8 £ 3 i , 3 3 3 8 j 8 5 8 8 E S 8 3 8 E 6 8 S 6 C 1 8 £ S l ( 8 6 8 8 8 8 £ 8 C £ 8 £ S 6 8 8 S 8 S 8 B 3 S S 8 8 8 8 8 3 

t 

$ 9 $| D9S53SSSS999933S9S9S9&S9SS99S9S9S99S99S9S9SS99S9399SS89S9SSSS9SSS93S9S999SS9 | 

1 : • < s 1 1 3 1 id 1: '.i •: 11 iii5 :i .1 n 11 r :: :i 1: :sz: zsiisasi a u :j 3 ! t;a<? •!!!:• ‘.it: •.< ::ar i>i: c; :: « 71 ists n u;)M 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*STRT 

1-5 

A5 


YMD 

11-20 

F10.5 

YYMMDD. of start time 

HMS 

21-30 

F10.5 

HHMMSS.SSS of start time 


4.2.7 


*1 r.OP Control Card 


The *7RQP control card is used to supply the program 
with executive . information for processing tropospheric 
refraction. If no card is present, each observation will 
be corrected for tropospheric refraction using a table of 
corrections based upon average met. data (see Section 5.1). 
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NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*TROP 

1-5 

A5 


ITROP 

6-10 

15 

refraction identifier 

0 = no refraction 

1 = refraction with 

met data supplied 

PRES 

11-20 

110.5 

pressure in mb. 

TEMP 

21-30 

110. 5 

temperature in °F 

VAPR 

31-40 

110.5 

vapor pressure in mb. 


ORIGINAL PAGE IS 
- OF POOR QUALITY 


4.2.8 *SKIP Control Card 

The *SKIP control card Is used to supply the program with the file 
number to be used on the RV tape (counting the tape label). If no card 
Is present, the second file is used. 


NAME 

COLUMNS 

FORMAT 

DESCRIPTION 

*SKIP 

1-5 

A5 


NFSKIP 

6-10 

15 

File number to be used for 
RV tape (counting the tape label) 


4.9 


2 September 1976 


4.2.9 *SIGH Control Card 

The *SIGH control card is used to supply the program with the value 
of CSIGH, the calculated standard deviation of the satellite height in 
meters. If no card is present, a value of 10 meters is assumed. 


*SIGH 1-5 A5 

CSIGH 11-20 F10.5 Calculated standard deviation of 

the satellite height 


4.10 


2 September 1976 
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S . 1 Tropospheric Refraction 

Each observation processed by the ARC computer program is 
corrected for tropospheric refraction although the program does 
not alter the value of the observation on the data tape. The 
program makes use of the Saastamoinen tropospheric refraction 
model as presented in the "Wallops Island Tropospheric Refraction 
Study and Analysis", C.C. Goad, 1974, except that a constant 
elevation angle of 90° is used. The correction is given by: 

dr - 0.002277 [P + (0. 05+1255. 0/T)V] 


where 

P is the surface pressure in mb. , 

T is the surface temperature in ®K, and 

V is the surface vapor pressure in mb. 

The ARC program contains a tabic of corrections based upon 
average values of P, T, and V for each month of the year and 
in latitude increments of 5° from -70° to +65°. This table 
is linearly interpolated at one second intervals to obtain 
the tropospheric refraction correction for observation 
residuals and sea surface height. The tropospheric refraction 
correction written and the output tape is evaluated at the 
time (latitude) of the smoothed observation. The interpolation 
tabic is given in Table f. . 1 . 
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5.2 


5.2 Altitude 


Each altitude observation processed by the ARC computer 
program is corrected for bias while computing residuals, 

SSHITE etc. The value of the bias correction is constant 
and is zero unless supplied by the user (see Section 4.2.1). 

The program accepts separate bias values for the global and 
intensive modes. 

5.3 Geoid 

During the computation of the altimeter observation 
residual it is necessary to calculate the modeled height of 
the geoid. The ARC computer program uses the Marsh-Vincent 
geoid (see "Geophysical Surveys", 1974, pp. 481-511, D. Reidel 
Pub. Co., Dordrecht, Holland). The geoid height is evaluated 
by linearly interpolating a table of values given in one-degree 
latitude by one-degree longitude increments. The interpolation 
is performed every second for observation residuals and sea 
surface height. The value written on the output tape is eval- 
uated at the time (latitude and longitude) of the smoothed 
observation . 

5.4 Tide 

During the computation of the altimeter observation 
residuals it is also necessary to calculate the modeled tide 
height. The ARC computer program uses the Hendershott tide 
model (see "Ocean Tides", M.C. Hendershott, EOS Trans., 

Febr. 1973) which is evaluated every second for residuals and 
sea surface heights. The value written on the output tape is 
evaluated at the time of the smoothed observation. 
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5.5 Ephemeris Interpolation 

In order to calculate observation residuals and sea 
surface heights, and evaluate the geoid and tide heights, it 
is necessary to interpolate the satellite latitude, longitude, 
and height above the reference ellipsoid. The ARC computer 
program uses a fixed mesh, fixed order Lagrange interpolation 
scheme, where the order is supplied by the user (see Section 
4.2.3) and can vary from 3 to 7 inclusive. The satellite 
latitude and longitude are evaluated every second. The values 
of latitude and longitude written on the output tape are 
evaluated at the time of the smoothed observation. The height 
above the reference ellipsoid is evaluated at every observation. 
The Lagrange interpolation scheme was chosen both because 
of its accuracy and simplicity (no difference tables are 
required and interpolation coefficients need only be computed 
once for a fixed mesh interpolation) . However it should be 
noted that the interpolation routine is not capable of 
extrapolating the epheineris data. Accordingly, the user should 
take care that several records of information are present 
on the ephemeris tape both before and after any possible 
observation (see Section 2.1). 


5.6 Orbit Geometry 

The geometry of the observation residuals, sea surface 
height, and associated parameters are given in Figure 5.1. 
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The computed observation residual is given by 
y = 0-C 

where y is the observation residual, 0 is the corrected 
observation, and C is the computed observation. The cor- 
rected observation is given by 
0 = O' - dr + B 

where O' is the actual observation, dr is the tropospheric 
refraction correction (see section 5.1), and B is the a priori 
altitude bias. The computed observation is given by 
C = H - G - T 

where H is the interpolated height above the reference 
ellipsoid (see section 5.5), G is the computed height of the 
geoid (see section 5.3), and T is the tide height (see section 
5.4) 

The computed sea surface height is given by 
SSH = H - 0 

and the computed smoothed sea surface height is given by 
SMOOSH = HS - SMOALT 

where SMOALT is the smoothed observation and HS is the height 
above the reference ellipsoid evaluated at the same time as 
SMOALT. 

The computed altimeter geoid is given by 
G’ = SSH - T. 
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5.7 Timing 

Each of t;._ parameters computed by ARC and written and 
the output tape (i.e., LAT, LONG, TREE, GHITE, THITE, SATHT, 
SSHITE, and SMOSSH) has a time of evaluation. The LAT, LONG, 
TREE, GHITE, THITE, and SMOSSH variables are evaluated at the 
time of the smoothed observation (at FODAY). If EODAY is not 
available from the GAP input tape, it is computed by ARC and 
LAT, LONG, and SMOSSH are set to -8888.0. SATHT, SSHITE, 

CALT, ALT, and the observation residuals computed by ARC are 
time tagged by the following algorithm: 

t i = FRAMTI + At + (i-l)fi 1 

where t^ is the time of the sample of SATHT, SSHITE, CALT, 

ALT, or the observation residual, At is the amount of time 
added to the frame time (FRAMTI) to obtain the time of the 
first sample, and 5t is the amount, of time elapsed between 
any two adjacent samples occurring in the same frame. 

The ARC program computes corrected values of At and 
fit for CALT and ALT and places them in the timing and header 

record. The corrections made are for transit time, filter 

( -I ^ 

delays and tracking loop delay, (see Tables 5.2, 5.3, and 5.4)'. 





PROGRAM MESSAGES 


ARC 

*** ARC NORMAL TERMINATION *** 

INIT 

DATA CARD NOT ACCEPTABLE - STOP 
HEADER 

GAP FLAG NOT SET 

PROCS 

*** NORMAL TERMINATION OF PROCESSING SECTION *** 
UNRECOGNIZABLE FORMAT - RECORD SKIPPED 
READER 

*** EOF READ ON OBSERVATION TAPE *** 

RE ADRV 

*** EOF READ GEODYN RV TAPE - STOP *** 

*** GEODYN RV TAPE CANNOT BE BACKSPACED - STOP ** 

TROP 

LATITUDE OF IS OUTSIDE TABLES 

GEOID 

SATELLITE LATITUDE = IS OUTSIDE OF GEOID - STOP 

PLTSAU 

DIMENSIONS IN PLTSAU INSUFFICIENT - STOP 

PLOTR 

NOTHING TO PLOT 
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Program Flowchart 



Print Output and ‘ 
■ Write Output Tape j 
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7.0 Product Jon Deck Setup 


The produ :t : on version of the ARC computer pro};r;irn can 
he executed by supplying the tape cards and the data control 
cards : 


$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


$ 

***EOF 


I DENT 112225, ARC 

OPTION . FORTRAN , RELMHM 
LOWLOAD 

ENTRY 

Binary decks and LINK cards. 
EXECUTE 

LIMITS 50 , 31 K , , 1 OK 

TAPE 0 1 , X 1 H , , nnnn , , RV 

TAPE 04 , X4P , , , , ARC 

TAPE 29.X9I), , minima , ,GAP 

FILE 0 3 , D 3 R , 2 0 S 

FILE 07.D7R.20S 

SYSOUT 1 1 

SYSOUT 12 

IJSF.RID PRODUCT IONS PRODUCTION 

PKMFL 02 , R , R, PRODUCTION/GEOI J) 

INCODH I BMP 

(Title Card) 

(Optional Control Cards) 

EX i) JOB 
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